IL NUOVO CIMENTO 



Vol. ?, N. ? 



On 
OS 

PL, 



On the parametric approximation in quantum optics 

G. M. D'Ariano, M. G. A. Paris and M. F. Sacchi 

Dipartimento di Fisica 'Alessandro Volta'and INFM — Unitd di PAVIA, 

Universita degli Studi di Pavia 

via A. Bassi 6, 1-27100 Pavia, Italy 

(ricevuto ?; approvato ?) 



> 

m 

o 

O 

o> 

ON 

a,: 

■*— > 

£ 



X 



Summary. — We perform the exact numerical diagonalization of the Hamiltonians 
that describe both degenerate and nondegenerate parametric amplifiers, by exploit- 
ing the conservation laws pertaining each device. We clarify the conditions under 
which the parametric approximation holds, showing that the most relevant require- 
ment is the coherence of the pump after the interaction, rather than its undepletion. 

PACS 03.65.-w, 42.50.-p, 42.65.Yj, 42.50.Ar . 



1. Introduction 

The interactions of different radiation modes through nonlinear crystals allow the 
generation of interesting states of light, which exhibit a rich variety of phenomena [1, 

2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]. Most of the theoretical analysis usually 
refers to situations where one mode — the so-called "pump" mode — is placed in a high- 
amplitude coherent state. In such cases the parametric approximation is widely used 
to compute the dynamical evolution [16]. In the parametric approximation the pump 
mode is classically treated as a c-number, thus neglecting both the depletion mechanism 
and quantum fluctuations. As a result, bilinear and trilinear Hamiltonians are reduced 
to linear and quadratic forms in the field operators, respectively, and hence some useful 
mathematical tools — typically decomposition formulas for Lie algebras — can be exploited 
for calculations [3, 17, 18]. 

In the validity regime of the parametric approximation different optical devices experi- 
mentally realize different fundamental unitary operators in quantum optics. For example, 
a beam splitter, by suitably mixing the signal state with a strong local oscillator at the 
same frequency, realizes the displacement operator, which generates coherent states from 
the vacuum. Similarly, a degenerate parametric amplifier realizes the squeezing opera- 
tor, which is the generator of the squeezed vacuum. Finally, a nondegenerate parametric 
amplifier realizes the two-mode squeezing operator, i.e. the generator of twin-beam. In 
the above scenario, it is matter of great interest to study the conditions under which 
the parametric approximation holds. This issue has been considered by a number of 
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authors [19, 20, 21, 22], however without giving a general validity criterion, which is the 
main concern of this paper. Quantum effects in two-mode optical amplifiers have been 
extensively analyzed [20, 21, 23, 24, 25, 26, 27, 28, 29, 30]. Frequency couplers with in- 
tensity dependent coupling have been studied [25, 26, 31], whereas the case of degenerate 
parametric amplifier has been considered by many authors [20, 21, 24, 27, 28, 29, 30]. 
Phase correlations [28] and the signal-pump degree of entanglement [29] have been ex- 
amined. The effect of pump squeezing has been also considered [30]. On the other hand, 
though trilinear processes have been thoroughly analyzed in Refs. [15, 24, 32], only little 
attention has been devoted to the parametric approximation in nondegenerate amplifiers 
[22]. 

The most explicit conditions for the validity of the parametric approximation have 
been derived in Refs. [19] and [20], for the beam splitter and the degenerate parametric 
amplifier, respectively. In both references sufficient conditions have been derived, which 
however can be widely breached in relevant cases of interest, as we will show in the 
following. 

In this paper we perform the exact numerical diagonalization of the full Hamilto- 
nians pertaining the three above-mentioned devices. As it was already noted by other 
authors [23, 24, 32] such a numerical treatment is made amenable by the presence of 
constants of motion that characterizes each Hamiltonian . In fact, the Hilbcrt space can 
be decomposed into the direct sum of subspaces that are invariant under the action of 
the unitary evolution. Therefore, one needs to diagonalize the Hamiltonian just inside 
each invariant subspace, thus considerably reducing the dimension of the diagonalization 
space [23, 24, 32]. 

We analyze in different sections the cases of the beam splitter, the degenerate para- 
metric amplifier and the nondegenerate parametric amplifier. The case of the beam 
splitter can be treated analytically, but we also present some numerical results in order 
to introduce the general approach that will be used for the parametric amplifiers. For 
each device, we look for the conditions under which the parametric approximation holds, 
for both vacuum and non- vacuum input signal states. The comparison between the exact 
evolution and the theoretical predictions from the parametric approximation is made in 
terms of the overlap O = y^Mj?outPth) between the state £> ut that exits the device 
and the theoretical state g^h obtained within the approximation. An explicit comparison 
in terms of photon number distributions and Wigner functions is also given for some 
interesting and representative cases. 

The main result of the paper is to show that the usual requirements for the validity 
of the parametric approximation, namely short interaction time and strong classical 
undepleted pump, are too strict. Indeed, we show that the only relevant requirement 
is the coherence of the pump after the interaction, rather than its undeplction. In fact, 
we will show typical examples in which the pump at the input is weak (one photon in 
average), after the interaction it is highly depleted, and notwithstanding the parametric 
approximation still holds. On the other hand, there are cases in which the pump after the 
interaction is only slightly depleted, however is no longer coherent, and the approximation 
fails. Finally, we show some interesting features such as pump squeezing and Schrodingcr- 
cat-like state generation that arise when the parametric approximation breaks down. 

2. Displacer 

The beam splitter is a passive device that couples two different modes of radiation at 
the same frequency through a first-order susceptibility-tensor \ medium. Such device 
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is widely used in quantum optics [1, 33], from homodyne/heterodyne detection [34], to 
directional couplers [35] and cavity QED [36]. In the rotating wave approximation and 
under phase-matching conditions, the beam splitter Hamiltonian writes in terms of the 
two mode operators a and b as follows 

(1) Hbs = k (arf + a)b) , 

where n is the coupling constant proportional to the x of the medium. The unitary 
evolution operator of the device in the interaction picture writes 

(2) U BS = cxp [-it (atf + a!b)] = e -' tanro6+ | cosr| a+a - 6t(, e- itanr ° t6 , 

where r is the interaction time rescaled by the coupling k. The factorization of the 
operator Ubs m Eq. (2) has been derived by applying the Baker-Campbell-Hausdorff 
formula for the SU(2) algebra [3, 17, 18]. The Heisenberg evolution of the field modes 
reads 

,„.. T y\ ( a \ - ( a cost — ib smr 

(3) U bs [ b ) U bs - y _ ia sinT + b C0ST 

From Eq. (3) it turns out that the transmissivity 9 at the beam splitter is given by the 
relation 8 = cos 2 r. 

The parametric approximation refers to situations in which one mode — say mode b — 
is excited in a strong coherent state. In this case in the first line of Eq. (2) the operator 
b might be replaced by a c-number, namely the complex amplitude (3 of the coherent 
state. Under this assumption, the evolution operator (2) would rewrite as the following 
displacement operator 

(4) D(-i0T) = exp [-it (/3a f + 0a)] . 

A more refined approximation that takes into account the 27r-periodicity in the exact 
Heisenberg equations (3) is 

(5) D(— i/3sinr) = cxp [— isinr (da/ + (3a)] . 

Indeed, the more precise result in Eq. (5) can be obtained by recasting the factorized 
expression in Eq. (2) in normal order with respect to mode a, after taking the expectation 
over mode b [37, 38]. The simple form of the bilinear Hamiltonian in Eq. (1) allows to 
clarify the conditions under which the parametric approximation (5) holds [19]. A set of 
sufficient requirements are given by 

|/3| — s- oo , sinr — ► 

(6) |/3| sinr = constant , 

without any assumption on the "signal" state for mode a. Hence, by combining a signal 
input state g- m with a strong coherent local oscillator \[3) in a beam splitter with very high 
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transmissivity, one can achieve the displacement operator in Eq. (5). The theoretically 
expected state g^ then writes 

(7) p^h = D(— 1/3 sin t) g- m D {— i/3sinr) . 

Here we present some numerical results concerning the exact unitary evolution of 
Eq. (2). The dynamics generated by the Hamiltonian (1) preserves the total number of 
photons involved in the process, in agreement with the following commutation relation 



(8) 



H BS ,a ] a + b ] b 



= 0. 



Therefore, it is convenient to decompose the Hilbert space TL a ® 7~Lb as a direct sum of 
subspaces with a fixed number N of photons, since these are invariant under the action 
of the unitary evolution operator (2). Such a decomposition can be written as follows 

(9) U a ®H h = (BJ^oHn 
where 

(10) H N = Span{|™}® \N - m) ,me [0,7V]} , 

Span{-} denoting the Hilbert subspacc linearly spanned by the orthogonal vectors within 
the brackets, and \n) ® \m) = \n, to) representing the common eigenvector of the number 
operator of the two modes. The decomposition in Eq. (9) makes the Hamiltonian (1) 
block-diagonal, namely 



(11) H BS = J2 h" . 



JV=0 



where Hn acts just inside the subspace TLm- Correspondingly, a generic two- mode state 
V>o) can be written in the orthogonal basis (10) as follows 



N 



(12) |V>q) = y^ y^ c m;A r_ m \m, N - to) 



jv=o m=0 



The diagonalization is performed inside each invariant subspace, and the truncation of the 
series in Eqs. (11) and (12) corresponds to fix the maximum eigenvalue of the constant 
of motion a^a + b^b. 

The state gout evaluated by the exact evolution operator (2) is given by 



(13) Qout = Tn[UBs(Qin®\P)(P\)U 1 



BSi 



where Tr?, denotes the partial trace on Tib. The comparison between the theoretical state 
g^h of Eq. (7) within the parametric approximation and the actual state g ut i s made 
in terms of the relative overlap 



(14) O = jTr[g th Q out 
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In the case of coherent input signal the overlap is evaluated analytically. One has 

(15) Ubs \ot) <8> \P) — |chcost — ifismr) |/3cosr — iasinr) , 
and thus 

(16) O = \{a — i/3sinr|acosT — i/3sinr)| = cxp ( — 4|a| 2 sin —J . 

From Eq. (16) it is apparent that the parametric approximation gives always exact 
results for vacuum input state (a = 0), whereas it is justified for coherent state as long 
as 4 1 a | 2 sin (r/2) <C 1, independently on the pump intensity. 

We introduce the quantity r* which represents, for a fixed value of the pump amplitude 
\P\, the maximum interaction time leading to an overlap larger than 99%. The value of 
r* clearly depends on the input signal: in agreement with Eq. (16) it is not defined for 
the vacuum (parametric approximation is exact), whereas for a coherent input signal |a) 
it is given by 

/ C 1 

t* = 2 arcsin A — - C = - (- In 0.99) 1/2 ~ 0.05 

VM 2 

The quantity r* also determines the maximum displacing amplitude \zm\ = |/3|sinr* 
that can be achieved by a beam splitter with a coherent pump \0). In Fig. 1 we have 
reported \zm\ for the vacuum, a coherent state and a number state as a function of the 
pump amplitude \[3\. The linear behavior of the plots indicates that r* is independent 
on the pump intensity. In the case of vacuum input we have complete energy transfer 
from the pump to the signal (slope of \zm\ y s \P\ equal to unit). Although they have the 
same energy, the coherent and number input states show different slopes, the coherent 
being more similar to the vacuum. Actually the set of coherent states is closed under the 
action of the displacement operator, so that the parametric approximation can fail only 
in predicting the exact amplitude of the output coherent state. 

We conclude that the first of requirements (6) is too tight. At least for coherent and 
number states, as long as the signal average photon number is less than the pump one, 
the beam splitter can "displace" the signal also for very weak pump. 

In the next sections we will deal with the problem of parametric approximation in 
nonlinear amplifiers. 

3. Squeezer 

The degenerate parametric amplifier couples a signal mode a at frequency cu a with a 
pump mode c at double frequency uj c = 2uj a . The interaction is mediated by the second- 
order susceptibility tensor x of the medium. Each photon in the pump mode produces a 
photon pair in the signal mode, giving rise to light with a number of interesting properties, 
such as phase-sensitive amplification, squeezing and antibunching [4, 5, 6, 7, 15]. In 
the rotating wave approximation and under phase-matching conditions the Hamiltonian 
writes 



(17) H DP = K(a 2 c 1 +a) 2 



c) 
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with k oc x ■ The corresponding unitary evolution operator in the interaction picture 
reads 

(18) U DP = exp [-it (a 2 c* + a* 2 c)] , 

where r represents a rescaled interaction time. The parametric approximation replaces 
the pump mode c by the complex amplitude (3 of the corresponding coherent state, such 
that the operator (18) rewrites as 

(19) S(-2i T (3) = exp [-it ((3a* 2 + Pa 2 )] , 

S'(C) being the squeezing operator [4]. In the case of coherent input signal |a), the 
predicted state at the output is the squeezed state 

(20) S{-2i T /3)\a) = S{-2iT(3)D(a)\0) = D(a)S(-2iT/3)\0) = \a , -2ir/3) , 

with a = acosh(— 2ir/3) + a sinh(— 2it(3). Notice that, differently from the beam splitter 
operator of Eq. (2), we have no method available to order Eq. (18) normally with respect 
to mode c [as in Eq. (2) for b] and then replace such mode by the c-number \/3\. Hence, 
we have no analogous nonpcrturbative method to estimate the validity of the parametric 
approximation. Hillery and Zubairy have been approached the question [20] in terms 
of a perturbation series for the propagator of the Hamiltonian (17). For initial vacuum 
state at mode a, they write the following conditions 

1/|/3| «1, r«l, 

(21) T e 4 ^ T « 1 , e 4 ^ T « |/3| . 

Here we evaluate the exact evolution generated by the operator (18) through numerical 
diagonalization of the Hamiltonian (17), using the method based on the constant of 
motion. In this case one has 



(22) 



Hdp, o o + 2c* c 



= 



and the Hilbcrt space TL a ®TL c is decomposed in terms of invariant subspaces correspond- 
ing to the eigenvalues of the constant of motion a) a + 2c' c, namely 

(23) Ti a ®Ti c = ®%™ H N , 
with 

(24) Ti N = Span {|7V - 2m) <g> \m) , m e [0, [N/2] }} , 

] denoting the integer part of x. Hence the Hamiltonian in Eq. (17) rewrites in the 
same fashion as in Eq. (11) and the block-diagonalization is performed for each hjsi , with 
N from to the maximum allowed value of the constant of motion. Similarly to Eq. 
(12), a generic two-mode state \iJ)q) is written as follows 

+ oo IN/2] 

(25) \ip ) = y^ ^ C N ^ 2m ,m \N - 2m, m) . 

N=0 m=0 
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The performances of a degenerate parametric amplifier in realizing the squeezing operator 
5(C) = cxp[l/2(Ca t2 - (a 2 )] of Eq. (19) are depicted in Fig. 2. In Fig. 2a we have 
reported the maximum interaction time t* that leads to an output signal whose overlap 
with the theoretical squeezed state is larger than 99%, as a function of the pump intensity 
|/3| 2 . In 2b we have shown the maximum squeezing parameter \Cm\ achievable by the 
amplifier, as a function of the pump amplitude \(3\. According to Eq. (19) one has 
Cm | = 2|/3|r*. In both pictures we have considered the vacuum, a coherent state and 
two different number states at the input of the amplifier. For the same set of input 
states, we have also shown in Fig. 3 the average signal photon number as a function of 
the interaction time r, for five different values of the pump amplitude. 

From Fig. 2 it turns out that the requirements for vacuum input signal in Eq. (21) 
are too strict. In particular, the two conditions in the second line are not satisfied for 
1 < \j3\ < 9 [see the line with triangles in Fig 2(a)]. Moreover, Fig. 2 shows that one can 
realize a squeezing operator even through a weak pump with just one photon. 

By definition the validity of the parametric approximation is guaranteed for t < t* , 
t* depending on |/3| and on the input state. However, we want to provide a general 
criterion that can be easily checked experimentally. As we will show in the following, the 
undeplction of the pump is not a valid criterion. We argue that the relevant parameter, 
in order to confirm whether the parametric approximation is justified or not, is the degree 
of coherence of the pump after the nonlinear interaction. These statements are supported 
by the following numerical results. 

Let us consider the case of a number input state \n = 1) with pump amplitude \/3\ = 9. 
From Fig. 2 one can extract the maximum interaction time t* ~ 0.073 for the validity 
of the parametric approximation, and the corresponding maximum squeezing parameter 
Cm | — 1-314. The average photon number of the output state can be drawn from Fig. 3c 
as (n)out — 9.68, corresponding to a pump depiction of about 5.4%. One might consider 
such a small depiction as the sign of the goodness of the parametric approximation. On 
the other hand, from Fig. 3c one recognizes the region 0.33 < r < 0.44, in which the 
output signal is even less excited than that in the above example, and consequently the 
pump is less depleted. Nevertheless the parametric approximation does not hold in such 
range of interaction time, since r is larger than r*. Let us now consider the Fano factor 
F = (An 2 )/(n) of the pump at the output. One finds that in the region 0.33 < r < 0.44 
the Fano factor is always larger than F — 1.13, whereas for r < r* it never exceeds 
F = 1.10. 

More generally, in all situations in which the parametric approximation is satisfied we 
found that the Fano factor of the pump at the output never exceeds F = 1.10. This holds 
also when the pump is weak (|/3| 2 = 1 -j- 10). Indeed, in this case the depletion of the 
pump can be strong, nevertheless the parametric approximation does not break down. 
In Fig. 4 we show the Fano factor of the output pump as a function of the interaction 
time r, for different values of the pump amplitude. Plots refer to vacuum input and to 
coherent state input \a = 1): similar plots can be obtained for other input states. 

As the condition of pump undeplction does not guarantee the validity of the paramet- 
ric approximation, so pump depletion by itself does not sign its failure: rather we have 
to consider the Fano factor of the pump. In order to stress this point, let us consider 
the extreme case of a pump with only one photon, and the input signal in the vacuum. 
The exact numerical solution indicates that the parametric approximation holds for in- 
teraction time up to t* ~ 0.42, the squeezing parameter and the output signal photon 
number increasing up to |Cm| — 0.84 and (n) ou t — 0.74, respectively [see Figs. 2 and 3a]. 
Correspondingly, the pump depletion grows up to 37% at r = t* . In spite of the strong 
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depletion, the pump preserves a good degree of coherence: the Fano factor achieves at 
most the value F = 1.10 at r = r*. 

In summary, the validity regime t < t* for the parametric approximation does not 
identify with the condition of pump undepletion, rather it corresponds to a Fano factor 
not exceeding the initial coherent level more than 10%. 

What happens beyond the parametric regime? For interaction time larger than t* 
new quantum effects arise at the output. In Fig. 5 we show the Wigner functions of both 
the signal and the pump modes at the output of the amplifier for t — t* and r = 2t*, 
with vacuum input and weak pump \(3\ = 1. In Fig. 6 the case of a stronger pump is 
given. As r increases, the pump first empties, then it starts refilling, preferably for even 
photon numbers, leading to oscillations in the photon number distribution. Remarkably, 
the corresponding Wigner functions of the pump and the signal exhibit interference in 
the phase space [39], the signal resembling a Schr6dinger-cat. 

4. Two-mode squeezer 

The nondegenerate parametric amplifier involves three different modes of the radiation 
field — say the signal a, the idler b and the pump c — which are coupled by a x^ nonlinear 
medium. The relation between the frequencies of the field modes is given by lu c = uj a +u)b. 
The Hamiltonian of the amplifier under phase-matched conditions can be written in the 
rotating wave approximation as follows 

(26) H NP = k (abc 1 + a^tfc) , 

with k oc x^ ■ The Hamiltonian in Eq. (26) describes also the case in which the frequen- 
cies pertaining modes a and b are the same, provided that the respective wave vectors 
and/or polarizations are different. The dynamics induced by the Hamiltonian (26) leads 
to a considerably rich variety of phenomena, such as generation of strongly correlated 
photon pairs by parametric downconversion [8, 9, 10], phase insensitive amplification 
[3, 16], generation of heterodyne eigenstates that are suitable for optimal phase detec- 
tion [11], polarization entanglement [12] and realization of Bell states [10, 12, 13, 14]. 
The unitary evolution operator in the interaction picture reads 

(27) U NP = cxp [-it (abc 1 + a^tfc)] , 

where r represents a rescaled interaction time. The parametric approximation replaces in 
Eq. (27) the pump mode c with the complex amplitude (3 of the corresponding coherent 
state, thus achieving the two-mode squeezing operator 

(28) S 2 {-iTJ3) = cxp [-it ({3rfb 1 + flab)] . 

The two-mode squeezing operator yields a suppression of the quantum fluctuations in 
one quadrature of the sum and difference of modes a ± b [40] . When applied to vacuum 
input, the unitary operator in Eq. (28) generates the so-called twin-beam 

DC 

(29) 5 2 (x)|0,0) = (l-|A| 2 ) 1 /2^ A "|n,n), 

n=0 
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where A = arg(%) tanh \x\- The expression in Eq. (29) can be easily derived by factorizing 
the §2 operator through the decomposition formulas for the SU(1,1) Lie algebra [3, 17, 
18]. 

The dynamics of the nondegenerate parametric amplifier admits two independent 
constants of motion. We choose them as follows 

(30) N= - [a t a + 6 t 6 + 2c t c] , K = a) a + c) c . 
Correspondingly, we decompose the Hilbert space Ha^Ti-b® H c in the direct sum 

(31) H a ® H b ® H c = ®^ =0 ©x=o Hnk , 
where the invariant subspaces "Hnk are given by 

Unk = Span {\K - to) ®\N - K - to) ®\m) , 

(32) to e [0, min(K, N - K)]} . 

The Hamiltonian Hnp and a generic three-mode state \ipo) will be consistently written 
as follows 

+ oo N 

(33) H NP = 2_^ 2_^ ^nk , 

N=0K=0 



+oc N min(K,N-K) 

(34) IV'o) = J^ J^ ^ CK-m,N-K-m,m\K -m,N - K -m,m) . 

N=0K=0 m=0 

To compute the exact dynamical evolution, one then diagonalizes each block Unk in 
Eq. (33) up to a fixed maximum value of N and makes the input state evolve in the 
representation of Eq. (34). 

As in the previous section, we have evaluated the maximum interaction time r* that 
provides an output state — in the signal and idler modes — whose overlap with the state 
predicted by the parametric approximation is larger than 99%. The time t* for vacuum 
input and number input \n = 1, n = 1} is plotted as a function of the pump intensity in 
Fig. 7a. The corresponding achievable two-mode squeezing parameter — the maximum 
argument \xm\ in the operator (29) — is represented in Fig. 7b. In Fig. 8 we show the 
average photon number (n) out of the output signal mode as a function of the interaction 
time, and for different values of the pump amplitude. Notice that the quantity a^a — b^b 
is conserved, so that, for the considered input states, the idler mode has the same average 
photon number as the signal one. 

As shown for the degenerate case, here also the requirement of a strong pump is not 
peremptory, whereas the undcplction of the pump does not guarantee the validity of 
the parametric approximation. Again, it is the Fano factor F of the pump after the 
interaction that well discriminates the working regimes of the amplifier. As long as F < 
1.10, the overlap between the states at the output and those predicted by the parametric 
approximation is larger than 99%. For interaction time longer than t* , the pump mode 
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reveals its quantum character, by showing oscillations in the number probability. This 
is illustrated in Fig. 9, where we report the photon number probabilities for both the 
signal and the pump modes in the case of vacuum input and pump amplitude equal to 
\/3\ = 5, for different values of the interaction time. 

5. Conclusions 

The quantum description of many optical devices is based on interaction Hamiltonians 
that couple different modes of radiation through the susceptibility tensor of the medium 
that supports the interaction. The theoretical predictions about such interactions are 
usually drawn in the so-called parametric approximation, i.e. by treating the pump mode 
classically as a fixed c-numbcr. Owing to such approximation, an analytical treatment 
is possible with the help of the factorization formulas for Lie algebras. 

In this paper we have investigated the conditions under which the parametric approx- 
imation holds in the treatment of x nonlinear amplifiers, by resorting to the exact 
diagonalization of the their full Hamiltonians. We have explicitly compared the states 
evaluated by the exact evolution with those predicted by the parametric approximation, 
in terms of the overlap between such states. On one hand, we have shown that the regime 
of validity of the parametric approximation is very large, including also the case of weak 
pump with 1 -r- 10 mean photon number. On the other, we have found that neither the 
condition of pump undepletion guarantees the goodness of the approximation, nor the 
condition of pump depiction signs its failure. We found that the degree of coherence of 
the pump after the interaction is a univocal parameter that discriminates the working 
regimes of the amplifiers. In terms of the pump Fano factor F we found that a deviation 
smaller than 10% guarantees an overlap larger than 99% between the states predicted 
within the parametric approximation and those evaluated by the exact Hamiltonian. 

For long interaction times the approximation breaks down, and the quantum character 
of the pump mode is revealed. Oscillations in the pump number probability appear and, 
correspondingly, the Wigner function of the signal mode assumes negative values and 
resembles a Schrodinger-cat state. 
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Fig. 1. - Performances of a beam splitter in achieving the displacement operator. We report the 
maximum displacing amplitude \zm\ achievable by a beam splitter as a function of the pump 
amplitude \/3\. In the picture triangles refers to vacuum input, squares to coherent state input 
| a = 1), and circles to number state input \n = 1). 
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Fig. 2. - Performances of a degenerate parametric amplifier in providing the squeezing operator 
5(C)- In both pictures triangles refers to vacuum input, squares to coherent state input \a = 1), 
circles to number state input \n = 1), and stars to number state input \n = 2). In (a) we 
report the quantity t* , namely the maximum interaction time that leads to an output signal 
whose overlap with the theoretical squeezed state is larger than 99%, as a function of the 
pump intensity |/3| 2 . In (b) we show the maximum squeezing parameter \£m\ achievable by the 
degenerate parametric amplifier, as a function of the pump amplitude \(3\. 
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Fig. 3. - Average photon number (n) ut of the signal at the output of a degenerate parametric 
amplifier, as a function of the interaction time r. In (a) the case of vacuum input, in (b) 
coherent input \a = 1), in (c) number input \n = 1), and in (d) number input \n = 2). Different 
line-styles refer to different pump amplitudes: j3 — 9 (dot-dot-dashed), j3 — 7 (dotted), j3 — 5 
(dot-dashed), (3 = 3 (dashed), /3 = 1 (solid). 
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Fig. 4. - Fano factor F of the pump at the output of a degenerate parametric amplifier, as a 
function of the interaction time r. In (a) the case of vacuum input; in (b) of coherent input 
\a = 1). Different line-styles refer to different pump amplitude: /3 = 9 (dot-dot-dashed), P — 7 
(dotted), p = 5 (dot-dashed), P = 3 (dashed), P = 1 (solid). 
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Fig. 5. - Contour plot Wigner functions of both the signal and the pump modes at the output 
of a degenerate parametric amplifier. The input signal is the vacuum, whereas the pump is 
initially in a coherent state /3 = —i. The time interaction is equal to r = r* = 0.42 in (a) and 
to r = 2r* = 0.84 in (b). 



16 



G. M. D'ARIANO, M. G. A. PARIS AND M. F. SACCHI 



Oh 



6 



C\2 

d 




"ft 



d 



C\2 

d 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

(c) - 




lllllllll. J, 



10 20 30 

n 



CQ 



I 

I 



C\2 



C\2 

I 

I 



"i — ' — r 



n — i — ' — r 
(b) 



J i L 




J i L 



-4-2 2 4 
Re a 




J i L 



J i L 



-4-2 2 4 
Re a 



Fig. 6. - Photon number probability and contour plot Wigner function for both the signal 
and the pump mode at the output of a degenerate parametric amplifier. The plots refer to a 
situation in which the signal mode is initially in the vacuum and the pump mode is excited 
to a coherent state with amplitude j3 — — 3i. In (a) and (b) the interaction time is equal to 
t = t* — 0.19, whereas (c) and (d) refer to an interaction time r = 0.43 > 2r*. In the first case 
the parametric approximation well describes the real interaction, which produces a squeezed 
vacuum state with squeezing parameter r — 1.146 corresponding to about 2 squeezing photons. 
On the other hand, parametric approximation does not hold in the second case. Notice that the 
break-down of parametric approximation is connected with the appearance of negative values 
in the Wigner function, which is a signature of quantum interference in the phase space. 
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Fig. 7. - Performances of a nondegenerate parametric amplifier in achieving the two-mode 
squeezing operator. In both pictures triangles refers to vacuum input and circles to photon 
number state |1, 1) input. In (a) we report the quantity r*, namely the maximum interaction 
time that leads to an output signal whose overlap with the theoretical state is larger than 99%, 
as a function of the pump intensity |/3| 2 . In (b) we show the corresponding maximum two-mode 
squeezing parameter \xm\ achievable by the nondegenerate parametric amplifier, as a function 
of the pump amplitude \/3\. 
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Fig. 8. - Average photon number (n)out of the signal at the output of a nondegenerate paramet- 
ric amplifier, as a function of the interaction time r. In (a) the case of vacuum input, and in (b) 
the case of photon number input |1, 1). Different line-styles refer to different pump amplitudes: 
/3 = 9 (dot-dot-dashed), (3 = 7 (dotted), (3 = 5 (dot-dashed), (3 = 3 (dashed), f3 = 1 (solid). 
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Fig. 9. - Photon number probabilities for both the signal and the pump modes at the output of 
a nondegenerate parametric amplifier. The plots refer to a situation in which the signal mode 
is initially in the (two-mode) vacuum and the pump mode is excited to a coherent state with 
amplitude j3 = — 5i. The interaction time is equal to r = t* = 0.214 in (a), and to r = 3r* 
in (b). In the first case the parametric approximation well describes the real interaction, which 
produces a twin-beam state with two-mode squeezing parameter \ = 0.789 corresponding to 
about 3.3 output photons. On the other hand, parametric approximation does not hold in the 
second case, as it can be easily recognized from the pump squeezing. The pump Fano factor in 
(b) is about F = 7.4. 



